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We discuss the statistical mechanics of magnetic flux lines in a finite-thickness slab 
of type-II superconductor. The long wavelength properties of a flux-line liquid in a slab 
geometry are described by a hydrodynamic free energy that incorporates the boundary 
conditions on the flux lines at the sample's surface as a surface contribution to the free 
energy. Bulk and surface weak disorder are modeled via Gaussian impurity potentials. This 
free energy is used to evaluate the two-dimensional structure factor of the flux-line tips at 
the sample surface. We find that surface interaction always dominates in determining the 
decay of translational correlations in the asymptotic long- wavelength limit. On the other 
hand, such large length scales have not been probed by the decoration experiments. Our 
results indicate that the translational correlations extracted from the analysis of the Bitter 
patterns are indeed representative of behavior of flux lines in the bulk. 
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1. Introduction 



The nature of the ordering of the magnetic flux array in the mixed state of high- 
temperature copper-oxide superconductors has received considerable experimental and 
theoretical attention in the last few years. It has been shown that fluctuations are im- 
portant in these materials and can lead to a number of new phases or regimes of the flux 
array, including entangled flux liquids, hexatic flux liquids [^0] and hexatic vortex glasses 
0,^. Most experiments probe the properties of the flux array indirectly by measuring 
bulk properties of the superconductors, such as transport, magnetization or mechanical 
dissipation. At present direct measurements of the microscopic order of the magnetic flux 
array are mainly limited to decoration experiments at low fields ■ These experiments 
aim to extract information on the vortex line configurations in the bulk of the material 
by imaging the pattern of the magnetic flux lines as they emerge at the surface of the 
sample. The surface patterns are determined by the interplay of thermal fluctuations and 
impurity disorder. Both these mechanisms can be responsible for disrupting translational 



and orientational order of vortex arrays [llO| , pT| . Surface roughness can also play a role in 



determining the surface magnetic patterns. It is clear that to interpret the experiments 
and assess whether one can indeed consider the surface patterns as representative of vortex 
line configurations in the bulk of the sample, one needs to understand what are the relative 
effects of bulk versus surface interactions and disorder in determining the configuration of 
the vortex tips as they emerge at the surface. Almost thirty years ago Pearl showed 



that the interaction between the tips of straight flux lines at a superconductor-vacuum in- 
terface decays as l/r± at large distances, with r± the distance between flux tips along the 
interface. In contrast, the interaction between flux-line elements in bulk decays exponen- 
tially at large distances. For this reason Huse has questioned the assumption that surface 
patterns are representative of flux lines configurations in the bulk and has argued that at 
low fields, where the intervortex separation is large compared to the penetration length, 
surface effects may play the dominant role in determining the magnetic fiux patterns seen 
at the surface [|TB| . 

By analyzing fiux decoration images one can extract quantitative information on the 
decay of both translational and orientational correlations of fiux line tips at the sample sur- 
face. In this paper we focus on the long-wavelength behavior of translational correlations 
in the fiux-line liquid phase. This case may be relevant to the interpretation of decoration 
experiments such as those by the AT&T group for the following reason. The decoration 
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experiments are carried out by quenching the sample in a smaU field from high to low 
temperature. The observed fiux patterns do not represent the equilibrium configurations 
of vortices at the low temperature where the decoration takes place, but equilibrium con- 
figurations corresponding to a higher temperature, Tf, where the fiux array falls out of 
equilibrium. While the value of Tf is not known, there are indication that it may be near 
the irreversibility line, Tim which in turn has been found to be very close to at low 
fields in the BSCCO samples used for the decorations 0. The experiments may then be 
probing a rather narrow range of temperatures between Tirr and Tf where the fiux array 



is in the polymer-like state proposed by Nelson |TJ] . The long wavelength static properties 
of such a glassy polymer can be described in terms of a hydrodynamic fiux-line liquid free 
energy [0. 

Translational correlations of fiux lines in three dimensions are described by the density- 
density correlation function. For convenience we consider the Fourier transform of this 
correlation function in the plane normal to the applied field. 



noS{q±, Zi,Z2) = < Sn{q±, zi) (5n(-q_L, Z2) > - < 5n{(i±,zi) >< 57i(-q^, 22) >, (1.1) 
where n(q_L,^) is the in-plane Fourier transform of the coarse-grained fiux-line density. 



N 

n{r^,z) = J2 5(rx-r,(z)). (1.2) 



Here ^^{z) is the position of the z-th vortex in the (x, y) plane as it wanders along the z 
(z II H) axis and dn{r±,z) = n{r±,z) — uq denotes the fiuctuation of the local density 
from its equilibrium value, no = B/(po, with B the average induction at equilibrium and 
(pQ the fiux quantum. The angular brackets denote a thermal average and the overbar the 
average over quenched impurity disorder. The subtracted term on the right hand side of 
Eq. ( p. . 1|) vanishes in the absence of quenched disorder. 

The two-dimensional structure function of a constant-z cross-section of fiux liquid is 
obtained from ([L1] ) by letting zi = Z2, 



noS2{q±,Zi) = < \Sn{q^±, zi)\^ > - \< 5n{q^±,Zi) >\ . (1.3) 

A factor of no has been extracted from the definition of the the structure function so that 
'S'2(g_L, 2i) ^ 1 as q'j, ^ 00, for all values of zi. 



The behavior of the fliix-hne structure function (LI) for bulk flux-hne hquids in the 
presence of both thermal fluctuations and quenched disorder has been discussed before fl^ . 
In an infinitely-thick superconductor the correlation of density fluctuations at different 
heights only depends on the distance, \zi — Z2\- The results are then most conveniently 
discussed in terms of the full three-dimensional structure function obtained by Fourier 
transforming ( pTTP with respect to zi — ^2, 



nQS{q^,q^) = < \5n{(\^, qz)\'^ > - I < '^?^(q_L,?z) > I • (1-4) 

For simplicity we only discuss in this section the contribution to the structure factor from 
thermal fluctuations, 5't(i?_l, Qz). This is easily calculated in the long wavelength limit 
from a simple hydrodynamic theory [P,|15|, with the result. 



ST{q±,qz) = ^ J ^ — 2 7 V' ^^•^> 

where CL{q±,qz) and 044(5^,^2) are the wave vector-dependent compressional and tilt 
moduli of a bulk flux- line liquid, as given for instance in Ref. In the long wavelength 
limit one can invert the transform by approximating the elastic constants in ( |1.5| ) with 
their values at = 0, with the result, 

SA<l^,z,-z,) = ^^f^e-l---l/^ii(-), (1.6) 
B3{q±)^\\{q±) 

where 



(1.7) 

Bs{q±) q± 

is the correlation length describing the decay of correlations in the z direction, with 
B3{q±) = CL{q±,0) and K{qj_) = C44{q±,0). The asymptotic long-wavelength {q± 0) 
behavior of the correlation function is determined by the constant values of the moduli at 
zero wavevector, -03(0) = ^ and K{0) = ^ -1-044(0), where 044(0) is the tilt coefficient of 



a single line at q^ = 0. [|T^. On the other hand, the nonlocality of the elastic constants is 
often important in the high-Tc superconductors even at small wave vectors because of the 
large values of the penetration lengths [|1^ . 



Correlations in any constant-^ cross section of the bulk are described by the two- 
dimensional structure factor defined in Eq. ( p. .31) and do not depend on zi. The contri- 
bution to this two-dimensional structure factor from thermal fluctuations is immediately 
obtained by letting zi = Z2 in Eq. (|1.6| ) , 



noksT 



UoKbT 

According to a well known sum rule that relates the value of the structure factor at zero- 
wavevector to the bulk modulus of the liquid, the denominator on the right hand side 
of the first equality of Eq. ( |1.8|) can be interpreted as an effective two-dimensional bulk 
modulus. This effective bulk modulus diverges in the long wavelength limit, due to the 
divergence of the correlation length ^y. As a consequence the two-dimensional structure 
factor of any constant-z cross-section of a bulk fiux-line liquid vanishes linearly as q±_ 0. 



As discussed elsewhere [|15| , p!8[ , this behavior arises here from the constraint that fiux lines 
cannot start nor terminate inside the sample. 

Decoration experiments measure correlations of flux lines in a constant-^; cross-section 
at the surface of a sample of flnite thickness L. Most researchers have implicitly assumed 
that these surface correlations are representative of correlations in a constant-^; cross- 



section of bulk. If this is the case Eq. ( |1.8|) should describes the correlations of flux-line 
tips at the surface as extracted from the decoration experiments. On the other hand, it is 
clear that whenever > L flnite-size effects may become important. Since {q±) diverges 
as Qj^ — > flnite-size effects may in fact dominate the long wavelength behavior. In addition 
we have mentioned that in a flnite size sample the boundary conditions modify the pair 
interaction between flux lines at the surface. The boundary conditions are determined 
by the requirement that on length scales larger than the penetration length the spatially 
uniform magnetic fleld outside the sample must equal the fleld of the vortex tips at the 
surface. In other words the superconductor surface can be thought of as the boundary 
between a strongly anisotropic magnetic medium where the fleld is concentrated in the flirx 
lines and an isotropic magnetic medium (the vacuum). At the surface the flux tips then 
behave like magnetic monopoles of "charge" (I)q/2tx and therefore interact at large distances 
via a repulsive Coulomb-like potential, 14 (rj^) = (^o/27r)^(l/rj^). The two-dimensional 
correlations of flux-line tips as they leave the sample could then be very different from those 
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in a constant-^ cross section deep in the bulk. If in fact the long-range surface interaction 



would dominate in determining long- wavelength surface properties ||13|| , then decorations 
would image the configurations of a two-dimensional liquid of vortices interacting via a 
l/r_L potential. The long- wavelength bulk modulus of such a two-dimensional liquid is 
-B2(?_l) — -B^/(47rq'j^). The corresponding two-dimensional structure factor is given by 
5'2t(<?_l) = nokBT/B2{q±) and again vanishes linearly as q± — > 0, but with a different slope 
from the result of Eq. (|1.8| ). A more detailed analysis is clearly needed to discriminate 
between these two possibilities. 

In this paper we describe a model that can be used to study the long wavelength 
properties of flux arrays in flnite-size samples incorporating the appropriate boundary con- 
ditions for the flux lines at the sample surface. We consider a slab of a uniaxial anisotropic 
type-II superconductor. The fleld is applied along the c axis of the superconductor which 
is chosen as the z direction. The slab is inflnite in the xy plane and has thickness L in the 
z direction. The starting point for the calculation presented here is a model hydrodynamic 
free energy that is a simplifled version of the full hydrodynamic free energy obtained in 
T^. In essence we neglect any nonlocality and spatial inhomogeneity in the z direction 



other than the presence of the sample boundaries. The free energy is then written as the 
sum of a surface free energy, that incorporates the boundary conditions, and the usual hy- 
drodynamic free energy of a bulk flux liquid. This free energy can be used as the starting 
point to study both how the coupling to the bulk affects the surface properties and how 
the bulk behavior is modifled by the presence of the boundaries. In this paper we focus on 
translational correlations at the surface. Our main result is the two dimensional structure 
factor at one of the surfaces of the sample. We flnd that the contribution to this structure 
factor from thermal fluctuations is given by, 

SM,.,L)= Jff (1.9) 

where B2^^ {q±: L) is an effective two-dimensional surface bulk modulus, 

B'2"{q^,L)=B2iq±)+Bsiq±)^\\{q±)F{L/^\\), (1.10) 
with B2{q±) ~ B"^ / {4TTq±) [^. The crossover function F is given by 

. _ i?3^||+i?2COth(L/e||) 

^^^/^ll)-5. + 53^„coth(L/ei,)- ^'-''^ 
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If the ratio K/B^ is large, as it can indeed be at the low fields used in the decoration 
experiments , one can distinguish three regimes in the behavior of the effective surface 
bulk modulus as a function of the in-plane wave vector gj^, 

The crossover wave vector q*^ is defined as the wave vector where the correlation length 
^11 equals the sample size, i\\{q\_) = L. For qj, — ^ the correlation length ^[| diverges and 
the fiux lines are essentially rigid. The surface interaction dominates in this limit and the 
surface bulk modulus is simply that of a 2d liquid of magnetic monopoles. When q^_> 
and > L, corresponding to qjL < q*^, fiux lines are still correlated in the z direction 
over the entire thickness of the sample, but the exponential interaction between straight 
fiux lines in bulk dominates in this case and the surface modulus is LB^. Finally, when 
^11 < L, 01 q± > q*^^ fiux lines are correlated over a length smaller than the sample 
size and the surface bulk modulus is determined by the compressional energy per unit 
area of an array of essentially straight fiux lines of length . The typical behavior of the 
surface structure factor as it crosses over from a linear function of q± at very small wave 
vectors to a linear function of q±_ with a smaller slope at large wave vectors is shown in 
Fig. 1. For both gj^ < 1/L and q±_ > q1_ [^, the surface structure factor ( |1.9|) is predicted 
to be linear in q±. The dependence of the slope on the areal density no of fiux lines is, 
however, different in the two regimes. For q± < 1/L we find S2t{(1±, L) ~ q±/no, while 

1 /2 

for q± > (j^, S2t{(1±tL) ~ q±/no ■ The latter behavior is consistent with experiments. 
The experiments are carried out on slabs which are typically of 1mm extent and 5 — 30/um 
thickness. Due to the finite size of the decoration images only information on the structure 
function for q± > 0.4/im can typically be extracted from the experiments that are therefore 
unable to probe the range q± < 1/L. In fact the experimental structure factors appear to 
level out in the small q± range [^. One can therefore conclude that the patterns probed 
by decorations are indeed representative of bulk behavior. Only the asymptotic behavior 
at very long length scales is dominated by surface effects. These large length scales have 
not, however, been accessible in experiments. 

In Section 2 we present a simple model hydrodynamic free energy for a fiux-line liquid 
in a finite-size superconductor sample that properly incorporates the boundary conditions 
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at the superconductor /vacuum surface. This free enrgy is used to evaluate the thermal 
part of the structure function at the sample's surface. In Section 3 we introduce disorder in 
the hydrodynamical treatment and discuss its effect on the translational correlations at the 
sample surface. The comparison with experiments is discussed in Section 4. In Appendix 
A we present the results for the full three-dimensional structure function defined in Eq. 



2. Hydrodynamic Free Energy and Surface Structure Function 

The pair interaction between flux lines in a finite-size superconductor sample has been 
derived before in the London approximation for both isotropic ||2^ and anisotropic [^|l 



superconductors. In Ref. [19 we calculated the magnetic energy of the flux-line array 
in a semi-inflnite sample of anisotropic Cu02 superconductor occupying the half space 
2; < 0, where the z direction is the c axis of the material and the fleld is applied along 
the c axis. The calculation of the energy for a superconducting slab of thickness L can be 
carried out along the same lines. In general the magnetic energy from the stray flelds in 
the region ouside the superconductor can always be rewritten as a surface contribution to 
the pair interaction between flux lines that decays exponentially with the distance from 
the interface. For straight lines the pair interaction behaves as discussed by Pearl and it 



decays as at large distances at the superconductor's surfaces. Following [|T9[, this 

pair interaction can then be used to obtain the coarse-grained hydrodynamic free energy. 
The presence of the superconductor/vacuum interface modifies the compressional and tilt 
elastic constants of the fiux-line array. In addition to the familiar nonlocality of the elastic 
constants associated with the range of the repulsive interaction between fiux-line elements, 
the presence of the interface introduces additional nonlocalities in the z direction. The 
corresponding surface contributions to the elastic constants depend exponentially on the 
distance of the deformed fiux volume from the interface. The surface contribution to the 
wave vector-dependent bulk modulus diverges as l/q±_ at small wave vectors, as expected 
for particles interacting via a potential in two dimensions. The surface contribution 
to the tilt modulus is negative and finite in the small wavevector limit. 

Instead of using this general hydrodynamic free energy, here we propose a much sim- 
pler model free energy that neglects all nonlocalities in the z direction, other than the 
presence of the superconductor /vacuum boundaries. For clarity we consider in this section 



the case of a clean superconductor. The effect of weak impurity disorder on the trans- 
lational correlation functions of a finite-thickness flux-line liquid will be discussed in the 
next section. The flux liquid free energy is then written as the sum of bulk and surface 
contributions, 

F = Fb + Fs, (2.1) 

where Fb is the usual hydrodynamic free energy for a bulk flux-line liquid. It includes 
terms quadratic in the hydrodynamic fields, which are the density field defined in ( |1.2| ) 
and a tilt field, defined as, 

= (2.2) 



dz 

1=1 



The bulk free energy is given by [|T] 



2n, 



^V/ d4i?3(g±)|<5n(q^,z)p + K(g^)|t(q^,z)|n. (2.3) 



Since we are allowing for in-plane nonlocality of the elastic constants, we have written the 
free energy in terms of the Fourier components of the hydrodynamic densities. There are 
two surface contributions to the free energy, 

= ^ E {52(g±)|5n,(qx, ^ = 0)|2 + S2(g±)|5n„(q^, z = L)|2}, (2.4) 



where B2{q±) is the surface contribution to the bulk modulus obtained in and given 
by 

^2(0^) = , (2.5) 

4rrg^ {I + qlXl,mq±Xa, + + qiK,) 

For q±Xab << 1 it is well approximated by the compressional modulus of a two-dimensional 

fluid of particles interacting via a potential, B2{q±) — B"^ /{4Tvq±). 

Statistical averages over vortex line conflgurations have to be carried out over the free 

energy ( |2.1| ) with the constraint that flux lines cannot start nor stop inside the medium, 

djn{ci±,z)+i(i±-t{(i±,z)=0. (2.6) 

It is convenient to separate t(qj^, z) in its longitudinal and transverse parts, t = q±tL+z x 
q^ty, with = q± • t and tr = {z x q±) -t. As a result of the constraint (|2.6|), the density 
and the longitudinal part of the tangent fleld are not independent hydrodynamic variables. 
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In fact they can both be expressed expressed in terms of a "vector potential" u(q±,z), 
which plays the role of a two-dimensional displacement field for the fiux-line liquid, 

5n(q_L, z) = -noiq± ■ u{q±, z) 

(2.7) 

^L(q_L, z) = no^^qx • u(q_L, z). 

The transverse degrees of freedom It are decoupled from the longitudinal ones. Since in 
this paper we are only interested in calculating correlation functions of density fluctuations, 
we only need to consider the longitudinal part of the free energy. This can be expressed 
entirely in terms of ML(q±, z) = c{±- u(q^, z) so that the constraint (|2.6| ) is automatically 
satisfied, with the result. 



^'"'-^T. I |^(^^)|^f + 53(,J,i|«.(qx,.)|^ 

° ^ (2.8) 

+ [d{z)+5{z-L)\B2{qi_)ql\uL{ci^,z)\''Y 

The tilt energy provides the coupling between bulk and surface terms. The model free 
energy of Eq. (|2.8|) has the same structure of model free energies used to study surface 



phase transitions, particularly in the context of wetting ||2^ . The longitudinal component 
of the vector potential takes the place of the order parameter and the presence of interfaces 
introduces spatial inhomogeneities in the order parameter. 

One could now proceed directly to evaluate correlation functions of the field by 
taking statistical averages with the free energy (|2.8|) . The two-dimensional structure factor 
of the flux-line tips at the top surface is simply given by S'2t(?_l) = ?^o?l < |^l(ci_L5 -^) P 
where the angular bracket denote the statistical average over the free energy ( ^.8| ). Since 
the free energy is quadratic in the fields, these correlation functions can be calculated 
exactly. On the other hand, due to the coupling of the field at the surface to the field in 
the bulk of the sample the calculation is somewhat lengthy and tedious. A much simpler 
way to obtain the same result is to use linear response theory. Let us apply a spatially 
inhomogenoeus surface pressure 5p(rj^) at z = L and consider the linear response of the 
system to this perturbation. The surface pressure couples to the density and the free 
energy of the corresponding perturbation is given by 



dFp = — [ dr±dp{r_i_)dn{rj_,L) 
no J 

= 5Z ^q^^Pi^i-hL{(l±, L). 



(2.9) 



9 



The response to this perturbation is a nonvanishing average displacement, < UL{q±, z) >p, 
that can be determined by minimizing the total free energy, that is, by requiring, 

f =0. (2.10) 

The minimization of the total free energy yields an equation for the displacement field in 
the bulk of the sample, 

-Kdj < ul{z) >p +qlBs < >p= 0, (2.11) 

for z 7^ 0, L, and boundary conditions for the displacement field at the superconduc- 
tor/vacuum surfaces, 

K[d^ < UL >p]^^L + ^i^2 < ul{L) >p -iq±6p = 0, (2.12) 
-K[d,<UL >p]^^^ + qlB2<UL{0) >p=0. (2.13) 
By solving Eq. ( p.ll|) with the boundary conditions (|2.12|) and ( |2.13| ) we obtain, 

< WL(q±, z) >p= Aie-^/^ii + Aae^/^ii , (2.14) 
where $.\\{q±) is the correlation length defined in Eq. ( |1.7| ) and 



' 2q^ [Sf^2^S2]sinh(L/^l|) + 2S3?l|S2COsh(L/^ll^ ^^"^^^ 



(2.15) 



' 2q^ [Sf + Si] sinh(L/^l|)+ 2^3^11 S2COsh(L/^l|) ^^"^^^^ 
The surface displacement is then given by, 

where F(L/^||) is the crossover function given in Eq. (|1 . 1 1|) . As expected, the response 
to the applied pressure is linear in the perturbation. The corresponding linear response 
function defined as 

5n{q±,L) -iq±no < UL{q±,L) >p 
XT{q±, L) = — — = — ^ , 2.17 

is the wave vector-dependent surface isothermal compressibility, which is in turn the re- 
ciprocal of the two-dimensional surface bulk modulus, B2^^ {q±, L) = l/x{q±, L). The 
denominator of Eq. ( |2.16| ) is then the effective surface bulk modulus given in Eq. ( |1.10|) . 
Finally, the surface structure factor is given by Eq. (|1.9|). 
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3. Weak Disorder 

The effect of weak disorder botli in the bulk of the sample and at the superconduc- 
tor/vacuum interfaces can be modeled in a standard way in terms of a random potential 
that couples to the density field. In general surface disorder may include surface roughness 
and therefore differ considerably from impurity disorder in bulk. We therefore model bulk 
and surface disorder separately by adding to the free energy given in Eq. ( |2.8| ) two terms, 

5Fd = Jdr±J^ {Vt{r^,z)5n{r^,z) + [5iz) + S{z-L)]Vs{r^)5n{r^,z)Y (3.1) 

where the random potentials Vfc(r^, z) and Vs{r±) represent the effect of random impurities 
and small scale inhomogeneities in the bulk and at the surface, respectively. If the defects 
are randomly distributed, as for instance in the case of oxygen vacancies, we expect that 
the quenched fluctuations in the impurity potentials will obey. 



H(r±, z)Vb{r'^, z') = A,5(ri - v'^)5{z - z'), 

(3.2) 

V;(ri)K(rY) = A,5(rx-rl), 



with Vb{Y^,z) = and Vs{y±_) = 0. 



It was shown by Nelson and Le Doussal [15] that in bulk flux-line liquids a hydrody- 



namic treatment of weak bulk disorder produces "Lorentzian squared" correction to the 
hydrodynamic result for the structure factor, given by 

(2 \ 2 

or, by inverting the Qz transform, 

SM{q±,z^ - Z2) = noAj „ , 7, , , ) [^\\{qi_) + \z^ - ^2|]e-l^^-^^l/«ii . (3.4) 

The contribution from weak disorder to two-dimensional density correlations in a constant- 
z cross-section of bulk, is then obtained from (^^) by letting zi = z^, 



, n|]Abg||(g^) 

Notice that this contribution to the two dimensional structure factor of a cross section of 
flux-line liquid also vanishes linearly with q±^ at small wave vectors. 
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To evaluate the two-dimensional structure factor of flux-line tips at the surface of 
a slab of type-II superconductor in a field, we proceed as in Section 2. The total free 
energy, F^^'> + Fd, of the flux- line liquid consists now of the sum of Eqs. ( p.8| ) and (|3.1|). 
To this we add the surface perturbation given in (^]^) and evaluate the response of the 
system by minimizing the total free energy. As in Section 2 we obtain an equation for the 
displacement field in the bulk of the sample, 

-Kdl < Ul{z) >p +qlB3 < ULiz) >p -iq±noVbici^, z) = 0, (3.6) 

for z 7^ 0, L, and boundary conditions for the displacement field at the superconduc- 
tor/vacuum surfaces, 

K[dz < ul >p ]^^^ + q\B2 < ul{L) >p -iq±noVs{q±) - iq±dp = 0, (3.7) 

-K[d^ < Ul >p]^^^ + qlB2 < ul{z = 0) >p -iq±nQVs{(\^) = 0. (3.8) 



The solution of Eq. (|3.6| ) with the boundary conditions ( p77|) and (|3.8| ) is given by 



<ULiq±,z) >p= A;e-^/?ii -fA'se^/^ii - j\z' cosh[{z - z') /^\\]Vb{ci^, z'), (3.9) 



with 



2q^ J + S|]sinh(L/^l|) + 2^3^11 i?2COsh(L/^ll)V^^^II ^2]5p(q±) 
+ noK,(q±) [(Ss^ii - B2) + (Sa^ii + S2)e^/«ii] 

+ noiBs^n - B2) I dzVb{(i±,z) cosh ^ ) -|- sinh f ^^-r- 

Jo L ^ ?|| ^ ^3^11 ^ 4l| 

= i + Bl] sinh(L/e|,)V 253^1,5, cosh(L/e,|) { t^^^" + ^^^'^^^^^ 



(3.10) 



+ noK(q±) [(Ss^ii + B2) + (Sa^ii - S2)e-^/«ii ] 
+ no(S3^ll + B2) [ dzVb{ci±, z) [cosh ^ ) + sinh f 



^11 ^ ^3^11 V 

(3.11) 

The solution (|3.9|) depends explicitly on the impurity potentials Vb and Vs and represents 
the longitudinal displacement field for a given realization of the disorder. In the presence 
of quenched disorder the fluctuation-dissipation theorem no longer holds and one cannot 



simply identify the ratio of the mean displacement < ^^(qj^,^) >p to the perturbation 
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Sp{q±) with the system's compressibihty. In fact the mean displacement in the presence 
of weak disorder is simply equal to that given in (|2.14|) for a clean superconductor. This is 
because the terms linear in the impurity potentials in Eq. ( |3.9| ) vanish when averaged over 
the quenched disorder and A[ = Ai and A'2 = A2. One then needs to calculate directly 
the correlation function of the fluctuations as defined in Eq. (|1.4|) , or, in terms of the 
linear response, < WL(q^,L) >p. 



S{q^, zi, Z2) = nog_L < UL{q±, zi) >p< ul{-c[±, Z2) >-< WL(q_L, zi) >p < UL(-q±, Z2) >^ 

(3.12) 

After some algebra, one finds that the two-dimensional structure factor at the surface 
2; = L of the slab, 5'2(g_L, L) = S{q±, z\ = L, Z2 = L), can be written as the sum of three 
contributions, 

5'2(g±, L) = 5'2T(g±, L) + S2dh{q±, L) + S'2sd(g_L, L), (3.13) 

where 5'2T(g±, L) is the contribution from thermal fluctuations that was already obtained 
in Section 2, 

S2T{q^,L)= Jff;^ . (3.14) 
Sf^(gi,L) 



with B2^^ {q±T L) the surface bulk modulus given in Eq. (|1.10|) . The other two terms in 



Eq. ( |3.13| ) represent the contributions from bulk and surface disorder, respectively. They 
are given by, 

S2bd{q±, L) — ^" 



2[i?f ^(g^,L)]^ [B2 + Ss^ll coth(L/e||)]' 

[Bl^l + El] coth(L/e||) + 252^3^11 + [BUI ' ^2] . , 

smh (L/,^11) 



(3.15) 



and 



C , n n^A. , i?|e^/sinh^(L/e||) 

{Bl^^{q^, L)f I {B2 + Sa^ii coth(L/^||)] ' 

While these expressions appear rather complicated, one can show that for all g^ of interest 
here the terms inside curly brackets in Eqs. (|3.15| ) and (|3.16| ) are of order one. 
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4. Discussion 



The surface structure factors extracted from the analysis of the low-field decoration 
experiments carried out by the AT&T group show clearly a large peak at q±_ = At: / (v^ao), 
where ao = (2/-\/3ao)^'^^ is the nearest neighbor vortex spacing in the triangular Abrikosov 
lattice, corresponding to the nearest- neighbor vortex spacing [Q]. At smaller wave vectors 
the structure functions obtained at different applied fields all decrease linearly with wave 
vector according to S2{q±-,L) ~ 0. At the smallest wave vectors probed in 

the experiments, the decay of the structure functions seem to be levelling out [|2^. This 
behavior is most apparent in the experiments at the lowest fields. 

To discuss the comparison of our results with these experimental findings, it is useful 
to consider some limiting forms of the expressions obtained in sections 2 and 3. 

When the sample thickness L is large compared to the correlation length the ex- 
pressions given in sections 2 and 3 for the surface structure factors reduce to those for 
the two-dimensional structure factors at the surface of a semi-infinite sample of type-II 
superconductor. In this limit, L » or q'j^ >> g^, the crossover function given in Eq. 
( imi) is F(L/^||) ^ 1 and we find 

S2T{q±, L) c:^ ^ q±, (4.1) 

In this limit all contributions to the surface structure factor are determined by an effective 



bulk modulus B2 + -Bs^y ~ Bs{l + ^/K/Bs)/q±, where we approximated B2 ~ Bs/q± 
T9[] . As indicated in Eqs. ( [4.1|) - (|4.3|), while both the contributions to the structure 



factor from thermal fluctuations and from bulk disorder are linear in q± in this regime, the 
contribution from surface disorder is proportional to q'j_. We argue below that the linear 
dependence of the surface structure factor on q± observed in the decoration experiments 



by the AT&T group 1^,^ probably corresponds to this regime. If our interpretation is 
correct, our results indicate that surface roughness does not play an important role in these 
experiments. 
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In the opposite limit, L « ^y, the flux hnes are essentiaUy straight over the en- 
tire thickness of the sample. The crossover function (|1 . 1 1|) becomes F{L/^\\) ^ {B2 + 
i?3L)/i?3^ll and the various contributions to the surface structure factor are given by, 

When — 0, i?2 ~ Q'J while S3 ~ (fj^ and. B2 always dominates in the denominator of 
Eqs. (ID-O, yielding, 

lim 5'2t(?±, L) ~ ^"^"^ ~ g_L, (4.7) 
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.0 ' 2So 



2 

2 



hm 52,,,(Q^,L)^li^~Qi. (4.9) 

If i?2 ~ B^/q^, the above limiting forms apply in the entire range q± « 1/L. Finally, if 
the ratio K/B^ is large enough so that there is a nonvanishing range of wavevectors where 
1/L < q'j^ < g^, in this range the surface structure factors level out to a constant value, 
given by, 

52r(g^,L)^^, (4.10) 



^2..(g±,L)^^, (4.11) 



^2 • 



S2sd{q^,L)^^j^. (4.12) 



We remark that if one evaluates the surface structure factors for a flnite-thickness super- 
conducting slab using free boundary conditions for the flux lines at the sample surfaces, 
one flnds that in the limit q± ^ Q the various contributions to S2{q±,L) have precisely 
the sample-size-dependent constant values given in ( ^.1C1| )-( ^.12| ). In fact for all values of 
q± where our hydrodynamic theory is relevant, the expressions for S2{q±,L) for the case 
of free boundary conditions on the flux lines are obtained from those given here by setting 
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i?2 = 0. The precise nature of the boundary conditions only affects the very smaU 
behavior of the correlation functions. 

We conclude by discussing in some detail the comparison of our findings with the 
experimental structure factors obtained by C.A. Murray and collaborators from low field 
decoration images [0,^ . The fields in these experiments are in the range 8 — 50G, corre- 
sponding to flux-line areal densities of 0.4 — 2.5 fxm~^. The smallest wave vectors at which 
one can construct a structure function from the decoration images is ~ OA^m-\ The 

1 /2 

long wavelength tail of 5*2 ((/x) is found to fit a q±/nQ in an intermediate range of wave 
vectors. The slope of S2{q±) in this linear region is typically ~ 5 x 10"'^ — 10~^^m. At 
the smallest wave vectors where data are available S2{q±) appears to be levelling out to a 
constant value ~ 0.02, almost independent of the applied field. 

Our model predicts a linear decrease of S2t{(1±) with qj_ in the limit q± 0. The 
corresponding slope is determined by B2 ~ ^0^0/ (^ttiJj,) according to Eq. ( |4.7D . One finds 
5'2t(?_l) ~ 2 X 10~^ fxm~^ {T /no)q±. Using T ~ 88K and no ~ l/xm~^, we obtain 
S2t{(1±) ~ 2 X 10~^/um q±, a slope over two order of magnitude smaller than oberved 
in experiments In addition the dependence of the slope on the density is not of the form 
obtained in the experimental fit. Finally, in the limit q± ^ both contributions to the 
structure functions from weak disorder vanish as q"^. All of this is consistent with the ex- 
pectation that this asymptotic long wavelength regime is never probed in the experiments. 



If the ratio a/ K/B^ < 1, then for all wave vectors accessible to experiments < L 
and the various contributions to the surface structure function are given by Eqs. (|4.1j )- 
( ^?3|) . Again we have a linear dependence of 5'2t(<?_l) on q±, but with a slope even smaller 
than in the asymptotic q± ^ case discussed above. A similar estimate can be carried 
out for the contribution from bulk disorder. In this case the size of the linear slope is 
determined by A;,. If the pinning is by isolated O2 vacancies, we can estimate ~ U^^^ab, 
with Uq ~ 5K a typical pinning energy barrier. The corresponding linear slope is then 
even smaller than for the thermal part of the correlation function. 

Our results are qualitatively consistent with the experimental findings if ^jKjBj, >> 
1. In this case we predict a crossover from the linear dependence of S'2(q'±) on q^ described 
by Eqs. ( ^?T| ) and ( [4.2| ) for qx_ > q]_ to the sample-size dependent constant values given 
by Eqs. ( [4.10| ) and ( [4.11] ). On the other hand, the predicted size of the constant value 
of S2{q±) for q± < and the linear slope at larger wave vectors are more than an 
order of magnitude smaller than from experiments, even when we take into account the 
in-plane nonlocality of the elastic constants. We can in turn try to fit our results to 
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experiments to extract experimental values for the elastic constants. For instance in Fig. 
1 by requiring S2{q± > 1/L) ~ nok bT / {BI'^'p L) ~ 0.02 for H = 8G, we find B^^'P/ks ~ 
70Kij,m~^ for L ~ 10/xm. This value is about 500 times smaller than -63(0) = B'^/4n, 
nor can the in-plane nonlocality of S3 account for the difference. We can then extract 
a value of the tilt modulus K by fitting the slope in the linear region, S2{q± > ~ 
{nokBT / ^K^^PBf^)q±_ ~ 0.02;um q^. We find K^'^P/ks ~ 4.4 x lO^kfim-^, a value not 
inconsistent with the theoretical value for K{0). The corresponding crossover wave vector is 
q]_ = (\ j V) ^ K^^P j B^"^ ~ l.nm~^. A mechanism that can account quantitatively for the 
very small value of S3 could be the strong downward renormalization of the compressional 
modulus that was predicted by Nelson and Le Doussal at low density ||T5[1 . 

In conclusion, our work indicates clearly that, while the decay of translational corre- 
lations at very long wavelengths {q± 0) is indeed governed by surface properties and 
finite size effects, the correlation functions extracted from the decoration experiments are 
in fact representative of bulk behavior. More data at small wave vectors and low fields 
with a detailed analysis of such data are, however, needed for a quantitative comparison. 

Work by M.C.M. was supported by the NSF through grant DMR-91-12330. D.R.N, 
acknowledges support from the NSF through the Harvard Materials Research Laboratory 
and through grant DMR-91-15491. We are grateful to C.A. Murray and D. Grier for 
sharing with me experimental results prior to publication. 

Appendix A. Three-dimensional Structure Function 

In this Appendix we discuss of the full three-dimensional structure function given in 
Eqs. ( p. . 1| ) . We only sketch the derivation of the thermal part of this correlation function, 
STiq±, z, zq). This derivation is instructive as an example of how to evaluate correlation 
functions directly by taking statistical averages with the free energy ( |2.(j| ) rather than by 
using linear response theory. 

It is convenient to expand the displacement field WL(q_L, z) in a Fourier series, 

WL(q_L,^) = -^|wo(q±) + 2^ [up{ci±)cos{pTTz/L) + Vp{ci±)sm{p7rz/L)]y (A.l) 

p=i 
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The free energy ( |2.8| ) is immediately rewritten in terms of the Fourier amphtudes of the 
displacement field, with the result, 



1 f °° 

2AL^^\ S «p(9±)(kp(q±)l' + l^p(q±)l') 



p= — oo 



p= — oo 



, ^2 2 



CXD 2 % 



p= — OO 



(A.2) 



with 



(A.3) 



where ?;o(<l±) = and we used that Up{—q±) = Up{q^±) and fp(— q±) = —Vp{q±). Using 
Eq. (|2.7| ), the thermal part of the three-dimensional structure function can be written in 
terms of the longitudinal displacement field. 



ST{q±, Zi,Z2) = nogi < ■UL(q±, 2:i)ML(-q±, ^2) > ■ 



(A.4) 



By inserting Eq. ( |A.1|) into Eq. ( |A.4|) we obtain 



2 00 00 

ST{q±, zi, Z2) = ^^^^^^ X X < Up{q±)u*p,{q_L) > cos{p7rzi/L)cos{p'7TZ2/L) 

p= — 00 p' = — oo 

+ < Vp{q±)v*,{q±) > sm{pTvzi/L) sin(p'7r2;2/L)|. 

(A.5) 

The correlation functions of the the Fourier amplitudes of the displacement field can be 
calculated by using standard tricks to deal with coupled gaussian integrals. The result is 

LkeT 



< Vp{q±)vp'{-ci_L) >= 5p^p' 



(A.6) 



and 



< Up{q±)Up>{-ci±) >=5. 



LkeT kBTB2ql 



ap{q±) apap, (1 + B2qlSoY - {B2qlS,Y 

X {(1 + B2qlS,)[l + (-1)^+^'] - B2qlS,[{-lY + (-1)^']}. 

(A.7) 

The quantities denoted by So and Se are special cases of more general sums that will be 
needed below. These are, 



X 1 cos{pTiz/L) 
^o[q±,z) = j 2^ 



1 cosh[(^ - L)/^|| 



ap{q±) Bs^\\ql sinh(L/^|| 
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with ^\\{q±) the correlation length given in Eq. (|1.7| ). The quantities So and 5*6 are given 
by the values of the above sums at 2; = 0, i.e., So = ^0(5^, 0) and S^ = Se(q'_L, 0). 
Finally, by inserting Eqs. (|A.71 ) and ( |A.6| ) into Eq. ( [A.5|) , one obtains. 



_ nokeT / cosh[(zi - Z2 - L)/g||] B2 

^^^^^^''^''>- S3^|| \ sinh(L/^||) B,+Bs^\i^i\F{L/^\^ 



X 



cosh[(2i + Z2- L)/^ii] ^ Bs^w cosh[(2:i - Z2)/^\\ 



sinh(L/e|| ) 52 + ^3^11 coth(L/e|| ) sinh^ (L/^n ) 

(A.IO) 

The two-dimensional structure factor in a constant-2; cross-section is obtained from Eq. 
( I^TOl) by letting zi = Z2, 



52^3^11 cosh[(22i-L)/e||] 



S2 + S3^l|Coth(L/^ll) sinh(L/^|| 



(A.ll) 



For zi = L and zi = this is identical to the espression given in Eq. ( |1.9| ). 

It is interesting to consider the correlation of density fluctuations on the two opposite 
surfaces of the sample, corresponding to the three dimensional structure function of ([A.lOj ) 
for zi = L and Z2 = 0. From Eq. (|A.10| ) we obtain 

^t(?±, 0) = 52t(?±, 0)i?(?i, L), (A.12) 

where S2t{<ILt^) = 5'2t(?±,-^) is the two-dimensional structure factor of one of the two 
surfaces, and 

^3^11 

^ S2sinh(L/^||)+S3^||Cosh(L/^||) ^^"^^^ 

measures the correlations between flux-line patterns at the two opposite ends of the sample. 
If the size of the sample is small compared to the correlation length, L « ^y, flux lines are 
straight throughout the sample and R{q±,L) ^ 1. The patterns on the two surfaces are 
perfectly correlated in this limit. Conversely, if L >> the patterns on the two surface 
are uncorrelated and R{q±, L) ^ 0. A measure of the deviation of R{q±, L) from 1 would 
give us information about the degree of flux-line wandering in the superconductor. 



19 



References 



[1] M.C. Marchetti and D.R. Nelson, Physica C 174, 40 (1991). 

[2] For a review see: D.R. Nelson, in Phenomenology and Applications of High- Temperature 
Superconductors, K.S. Bedell et al. eds. (Addison- Wesley, 1992). 

[3] D.S. Fisher, M.RA. Fisher, and D. Huse, Phys. Rev. B43, 130 (1991). 

[4] D.S. Fisher, in Phenomenology and Applications of High- Temperature Superconduc- 
tors, K.S. BedeU et al. eds. (Addison- Wesley, 1992). 

[5] P.L. Gammel, D.J. Bishop, G.J. Dolan, J.R. Kwo, C.A. Murray, L.F. Schneemeyer, 
and J.V. Waszczak, Phys. Rev. Lett. 59, 2592 (1987). 

[6] G.J. Dolan, G.V. Chandrasekar, T.R. Dinger, C. Field, and F. Holtzberg, Phys. Rev. 
Lett. 62, 827 (1989). 

[7] C.A. Murray, P.L. Gammel, D.J. Bishop, D.B. Mitzi, and A. Kapitulnik, Phys. Rev. 
Lett. 64, 2312 (1990). 

[8] D.G. Grier, C.A. Murray, C.A. BoUe, P.L. Gammel, D.J. Bishop, D.B. Mitzi, and A. 
Kapitulnik, Phys. Rev. Lett. 66, 2270 (1991). 

[9] P.L. Gammel, in Phenomenology and Applications of High- Temperature Superconduc- 
tors, K.S. Bedell et al. eds. (Addison- Wesley, 1992), and references therein. 
[10] M.C. Marchetti and D. R. Nelson, Phys. Rev. B41, 1910 (1990). 
[11] E.M. Chudnovsky, Phys. Rev. B43, 7831 (1991). 
[12] J. Pearl, J. App. Phys. 37, 4139 (1966). 
[13] D. A. Huse, Phys. Rev. B46, 8621 (1992). 
[14] D.R. Nelson and S. Seung, Phys. Rev. B39, 9153 (1989). 
[15] D. R. Nelson and P. LeDoussal, Phys. Rev. B42, 10113 (1990). 

[16] As pointed out in Ref [§, the single-line tilt coefficient contains a term that survives 
even in the limit of very large anisotropy and was missed in much previous literature. 
Using the results of [g, one obtains 044(0) = {nQ(j)l/167r'^Xli^)[ln{-fK)/-f'^ - {l/{2-f'^)) + 
(1/2)] . Here A^f, and Ac = ■jXab are the penetration lengths in the ab plane and along 
the c axis, respectively, with 7 the anisotropy parameter. Also, k = Xab/Cab, with ^ab 
the coherence length in the ab plane. 

[17] The reader may question the consistency of approximating the elastic constants by 
their values at = 0, while retaining their dependence on q±. It turns out that the 
main effect of the nonlocality on the properties of interest here is that of reducing the 
value of the bulk modulus at q± ~ l/ao? as compared to its value at q± = 0. This 
effect is retained in the approximation used here. 

[18] D.R. Nelson, Physica A177, 220 (1991). 

[19] M.C. Marchetti, Physica C 200, 155 (1992). 



20 



[20] The highest fields probed in the decorations are of the order of lOOG. For apphed 
fields along the c axis diamagnetization corrections are large in the thin samples used 
in the experiments. As a result S ~ i?, even for H < Hd- 

[21] We note that since we are using a hydrodynamic theory we are always discussing wave 
vectors q± < Go, where Go = 27r(2/-\/3no)^/^ is the shortest reciprocal lattice vector 
of the triangular Abrikosov lattice, corresponding to the location of the first peak of 
the two-dimensional structure function. 

[22] C.A. Murray, private communication. 

[23] E.H. Brandt, J. Low Temp. Phys. 42, 557 (1981). 

[24] O. Buisson, G. Carneiro and M. Doria, Physica C 185-189, 1465 (1991). 
[25] R. Lipowsky, Phys. Rev. Lett. 49, 1575 (1982). 



21 



Figure Captions 

Fig. 1. The thermal contribution to the surface structure function, as given in Eq. ( |1.8|) , 
for H = 8G (for this value of the field the first maximum of 5*2 (qj.) is at q± ~ 
4.26/um~^), T = 88K, Xab = 0.3/um, 7 = 55 and L = 25fxm. Here -03(0) was 
treated as a parameter determined by fitting our results to the small wave vector 
part of the data of C.A. Murray et al. as discussed in Section 3. 



22 



